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Effects caused by the electron electric dipole moment (EDM) in gadolinium garnets are considered. 
Experimental studies of these effects could improve current upper limit on the electron EDM by 
several orders of magnitude. We suggest a consistent theoretical model and perform calculations 
of observable effects in gadolinium gallium garnet and gadolinium iron garnet. Our calculation 
accounts for both direct and exchange diagrams. 



. I. INTRODUCTION 

o : 

Violation of the combined symmetry of charge conjugation (C) and parity (P) has been discovered in the decay 
of the K meson about 40 years ago |l| . The exact origin of this symmetry violation remains an enigma, although 
the so-called standard model of electroweak interactions can describe these processes phenomenologically. It has also 
been proposed by Sakharov that the matter-antimatter asymmetry observed in our Universe could have arisen 
from a CP-violating interaction active at an early stage of the Big Bang. The CP-violation implies a time-reversal (T) 
asymmetry, because there are strong reasons to believe that the combined CPT-symmetry should not be violated • 
Violation of the time reversal symmetry has been observed recently in decays of K-mesons Q . T- violation together 
with well known parity (P) violation (T,P-violation) provides a nonzero electric dipole moment (EDM) of a system 
in a stationary quantum state. This is why searches for EDM of elementary particles, atoms and molecules are very 
• important for studies of violations of fundamental symmetries . In the present work we concentrate on the electron 
EDM. The present best limitation on the electron EDM comes from the experiment with an atomic Thallium beam 

^: @, 

O 1 4 < 1.6 x 10" 27 e ■ cm. (1) 

m ; 

There have been recent suggestions 0, H| to improve the sensitivity to the electron EDM substantially by working 
with solids. In particular gadolinium gallium garnet, GdsGasO^ (GGG), and gadolinium iron garnet Gd3FesOi2 
(GdIG) have been suggested. The idea for searches of the electron EDM in solid state experiments was first suggested 
by Shapiro in 1968 |9j. There are two ways to perform such an experiment. The first one is to apply an external 
magnetic field to the solid. Polarization of electrons by the magnetic field causes alignment of electron EDMs, and 
hence induces a voltage across the sample that could be detected. Another possibility is to apply a strong external 
electric field. This would align the EDMs of bound electrons and hence lead to a simultaneous alignment of the 
' electron spins; the magnetic field arising from this alignment could be detected experimentally. An experiment of 
, this kind has been performed with nickel-zinc ferrite [lOj. Due to experimental limitations the result was not very 
impressive. However, according to estimates presented in Refs. 0,01 application of novel experimental techniques 
and measurements with GGG and GdIG makes this direction highly promising. 

The first calculation of the T,P-odd effects in GGG and GdIG has been performed in Ref. ( see a l so Refs. 
|l'2l IT^I. The mechanism considered in that work was similar to the mechanism of T,P-violation in atoms due to 
the T,P-odd nuclear forces [T^. In essence the Gd 3+ ion has been treated as a large nucleus that has some effective 
Schiff moment. External electrons that belong to 2 ~ ions penetrate inside the Gd 3+ ion and interact with the 
"Schiff moment" of the ion. This results in an effective interaction between the lattice deformation and the electron 
EDM. It has been also pointed out in Ref. that there is a contribution that can not be reduced to the "Schiff 
moment mechanism" , this contribution is due to the exchange between external and internal electrons. However 
the exchange diagrams estimated in Ref. pd| gave a small contribution. In the present work we have found a new 
class of exchange diagrams that are very important. Calculation of the effect with account of exchange diagrams is 
not a simple problem because in this case the logic that leads to the Schiff moment 0, IT3 is not valid. For the 
same reason the simplistic model of the electronic structure of Gadolonium garnet used in is not sufficient for 
the present calculation. Therefore to perform the present calculation we have developed a more accurate model that 
treats gadolinium and o xyg en electrons simultaneously. We believe that the present calculation is more accurate than 
that performed in Ref. 11]. Nevertheless, surprisingly, the results are very close. 
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II. THEORETICAL MODEL FOR ELECTRONIC STRUCTURE 



A. Single-electron effective potential 



Compounds of our interest, gadolinium iron garnet and gadolinium gallium garnet, are ionic crystals, consisting 
of Gd 3+ , Fe 3+ (or Ga 3+ ) and O 2- ions. Uncompensated electron spins are localized at Gd 3+ ion, which has a 4/ 7 
configuration, and Fe 3+ ion, whose electronic configuration is 4d 5 . Both paramagnetic ions contribute to the T,P-odd 
effect. However, contribution of the ion to the T,P-odd effect scales as Z 3 , where Z is the nuclear charge (see, e.g. 
Therefore, we neglect the contribution of Fe 3+ and consider only Gd 3+ ions. 

To describe an isolated Gd 3+ ion we use the effective potential in the following parametric form 

r (f + 7 1 r) 2 (e^r- + 1) r 

Here we use atomic units, Z is charge of the nucleus, Zi is charge of the core of the ion, and fi, d and -q are some 
parameters that describe the core. Solution of the Dirac equation with the potential J3J) gives wave functions and 
energies of electron states. We use the following values of the parameters 

^ = f.OO, d=1.00 77 = 2.35. (3) 

These values provide a fit of experimental energy levels for Gd 3+ (Zi — 4) and Gd 2+ (Zi — 3). Calculated and 
experimental [15J energy levels averaged over fine structure are shown in Table [I] 



Ion 


electron state Calculaton Experiment 


Gd 3 + 


4/ 


-363 


-355 




5d 


-281 




Gd 2 + 


5d 


-156 


-157 




6s 


-120 


-121 




6p 


-158 


-156 



TABLE I: Calculated and experimental energy levels with respect to the ionization limit. The levels are averaged over fine 
structure. Units 10 3 cm -1 . 



B. Account of Gd 3+ environment 



In the garnet structure each gadolinium ion is surrounded by eight oxygen ions O 2 in a dodecahedron configuration 
(resembles distorted cube), see Ref. 16]. The configuration at two different angles of view is shown in Fig. ^ The 




FIG. 1: Dodecahedron configuration of O 2 ions around Gd 3+ ion in garnet structure. 

Gd— O distance is r Q = 4.53ets, where as is the Bohr radius. We need to know wave functions of oxygen electrons 
inside the Gd 3+ ion. Electronic configuration of 2 ~ is I s 2 2s 2 2p 6 . Consider 2p orbitals of the ion. It can be shown 
that orbitals do not contribute to the effect and we need to consider only 2p a orbitals, pointing towards Gd 3+ . 
Every 2p a orbital is double occupied, so we have 16 electrons in the vicinity of Gd 3+ . In first approximation one can 
neglect distortion of the oxygen cube and use a linear combination of oxygen 2p a orbitals with definite symmetry with 
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respect to the cubic group. We are particularly interested in the \S) wavefunction, which is symmetric with respect 
to reflection of cube's main axes, and \P X ), \P y ) and \P Z ) wavefunctions, which change sign with reflection of x, y and 
z axes, correspondingly. They are of the form: 



\b) 


|1) 


+ |2) 


+ |3>- 


- K) + |5) + 


|6) 


+ |7) + |8) 








V8 






\P*) 


_ |1) 


+ |2) 


-|3>- 


- |4) + |5) + 


|6) 


- 17) - |8) 








V8 






\Py) 


_ |1> 


-|2) 


-|3>- 


- |4> + |5) - 


|6) 


- I?) + |8) 








V8 






\P.) 


_ |1> 


+ |2) 


+ |3>- 


h |4) - |5) - 


|6) 


- 17) - |8) 








V8 







(4) 

Here \n) denotes 2p a orbital of n-th oxygen ion as they are enumerated in Fig. El There are also three D-wave states 
and one F-wave state combined from oxygen 2p a orbitals, see Ref. but these states do not contribute to the 
effect and we neglect them. In the work the wave functions (@J have been matched to the gadolinium 6s and 6p 
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FIG. 2: Gd ion with surrounding eight oxygen ions. 



wave functions at the matching sphere of radius R ~ 2.5as around Gd. This is a simple and reliable way to describe 
penetration of oxygen electrons inside Gd 3+ . However, unfortunately this is not a consistent quantum mechanical 
description because the space is divided in two parts — inside and outside of the matching sphere. Our analysis 
has shown that calculation of exchange diagrams in this simple model is ambiguous. The point is that the effect 
arises in the third order of perturbation theory (see below). In such a high order one must use a consistent quantum 
mechanical description. The best solution of the problem would be a 3-dimcnsional Hartree-Fock or 3-dimensional 
pseudo-potential calculation for the cluster consisting of Gd and eight oxygen ions (Gd0 8 -cluster). However this is 
a very complex calculation. In the present work we suggest an intermediate solution, essentially a jelly model. We 
assume that eight oxygen ions provide a spherically symmetric attractive potential for electrons 



Vo{r) = -A e 



(5) 



where r a = 4.53 is Gd-0 distance, and A Q and D are some parameters. So effectively we replace the dodecahedron 
in Fig. ^ or the cube in Fig. |3 by a sphere. Thus the electrons are moving in the combined spherically symmetric 
potential 



V(r) =V Gd {r) + V {r), 



(6) 



where Vcd is the potential of Gd 3+ core (j2J , and Vo is the combined potential of oxygens JHJ . Solution of the Dirac 
equation with potential © gives single particle orbitals. In this picture the electronic configuration of the GdO 8 - 
cluster is ls 2 ...5s 2 5p 6 4/ 7 6s 2 6p 6 . Out of this ls 2 ...5s 2 5p 6 4/ 7 are the Gd 3+ electrons and 6s 2 6p 6 are in essence the 
oxygen electrons. So we have eight oxygen electrons, this is exactly what one needs to describe S and P 2p a states 
l|4|l. In this model one cannot describe D and F states combined from 2p a oxygen orbitals, but fortunately there is 
no need in these states because they do not contribute to the effect. 

How to determine the constants A a and D in the effective oxygen "jelly" potential (JSJ? The size of the potential D 
is not very important, it is clear that D ~ 1 (we remind that we use atomic units), and we set D = 1. The depth of 
the potential A Q is crucially important. Wave functions of oxygen ion 2 ~ have been calculated in Ref. 0, so ip2p„ 
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is known. We chose the depth A Q from the requirement that the wave functions ip6s and i/jg p satisfy the following 
conditions 

\As(R)\ = \^ 2p Aro-R,cos6 = l)\ , 

|V> 6 pCR,cos0= l/y/3) | = | V>2 P(r (r cos 0= 1)1 , (7) 

where R w 2.5. This is just an alternative formulation of the idea of dual description at r sa R, see Refs. [Til Il7j. 
There are two conditions J7J and only one parameter A a , so strictly speaking one cannot satisfy both conditions 
exactly. However at A a = 0.9 each of the conditions (JZJ) is satisfied with accuracy ~ 15%, and this is the value of A a 
which we use in our calculations. Thus, parameters of the "oxygen" potential (JjJJ are 

r Q = 4.53, A = 0.9, D = 1. (8) 

Let f(r) be the upper component of the Dirac spinor, ip(r) = f(r)/r ■ Yi m (9, </>). The electron density |/(t")| 2 for 6s, 
6f>i/2) an d 6P3/2 states calculated in the potential @ is shown in Fig. The density is peaked at r = r a = 4.53, 




FIG. 3: Radial electron density \f(r)\ 2 , V'fV) = f(r)/r-Yi m (8, 0), versus distance from the Gd nucleus. The solid line corresponds 
to the Qs-state, the dashed line corresponds to the Qp 1 / 2 -state, and the dotted line corresponds to the Gp^^-state. 

this is the effective description of oxygen ions. Energies of 4f, 5d, 6s, and 6p states in the potential © are listed 
in Table ITU The energies are given with respect to the ionization limit. It is worth to note that according to this 



State 


4/s/2 


4/7/2 


5d 3 / 2 5d 5/2 


6s 6pi/ 2 6p 3 / 2 


Energy 


-366 


-360 


-299 -295 


-311 -289 -286 



TABLE II: Energy levels in the potential © with respect to the ionization limit. Units 10 3 cm 1 . 

calculation the 5d-4f splitting is AEf d ~ 66 x 10 3 cto _1 . This is different from the splitting for an isolated Gd 3+ ion, 
AEfd w 100 x 10 3 cra -1 , see Ref. [ljj- However, the value AEfd ~ 66 x 10 3 cm _1 agrees well with that extrapolated 
for the Gd 3+ ion in the garnet environment, see Ref. [Isj . 



III. T,P-ODD ENERGY CORRECTION RELATED TO THE LATTICE DEFORMATION 

Similar to the approach 1 1 1 1 ] in this section we consider an external deformation of the lattice consisting in a shift 
of the Gd 3+ ion with respect to surrounding oxygen ions. Later we will relate the deformation to observable effects. 

We have already mentioned the T,P-odd effect arises in the third order of perturbation theory. So, there are three 
perturbation theory operators we have to consider. First of all this is the operator of the T,P-odd interaction of the 
electron EDM d e with atomic electric field E, see e.g. Ref. 0, 

V d = E . (9) 

Here 7 Q and X = 7o757 are Dirac 7-matrices. Because of the SchifPs theorem |l9f it is crucial to account for complex 
many-body screening effects, when working with the Hamiltonian 10 • Technically this means that the many-body 
perturbation theory with operator @ is very poorly convergent. The standard way [j| to avoid this complication is to 
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split the Hamiltonian © into two terms: Vd = — <i e 7 X • E = — d e S • £ — d e (j — 1)S • .E. Then due to the Schiff's 
theorem contribution of the first term to an observable effect is identically zero, so one can reduce the interaction 

V d -> VJ = -dei^o - 1)S • E. (10) 

The perturbation theory expansion with this operator is reasonably convergent. 

The second perturbation operator is related to the shift of the Gd 3+ ion with respect to the surrounding oxygen 
ions. Let us denote value of the shift by x. In our model, this corresponds to the shift of spherically symmetric Vo{r) 
(J5J along one of the axes, say z-axis, see Fig. 0] 



A Z A Z 




FIG. 4: Schematic picture, illustrating the shift of Vo(r) due to the lattice deformation. 

Therefore 

Voir) -> Vo(r) = V (r + x) = V (r) + {x ' r) = V (r)+xcos9^ . (11) 

r or or 

We keep the term linear in x, is the azimuthal angle. Thus the perturbation operator related to the lattice 
deformation reads 

V x (r) = x cos 9^ = -2x cos 6 ^ ~ fo) V Q (r) (12) 
Or D z 

The third perturbation operator is residual electron-electron Coulomb interaction, which is not included in the 
effective potential 

^,rj) ^ = E ^^SnWUTi) (13) 

Here and rj are radius-vectors of interacting electrons. 

Formula for the energy correction in the third order of perturbation theory reads, see e.g. Ref. [2(j 

E n '=2^ 2^ 75- ~ ^ 2^ 75-2-. ( 14 ) 

m k W mn Wkn m "> U nm 

where 

V = VJ + V x + V c . (15) 

In Eq. I|14f> we need to consider only the terms that are linear in each of the operators VJ (|10fl . V x l|12(l . and Vc i|13[l . 
There are 15 diagrams corresponding to Eq. I|14|l . these diagrams are presented in Fig. [SJ Each diagram contributes 
with a coefficient shown before the diagram. Summation over all intermediate states \k) and \m) and over all filled 
states \n) is assumed. The first four are the direct diagrams that correspond to the mechanism considered in |ll| . 
All other diagrams are exchange ones and therefore they contribute with sign (-). This is not so with respect to the 
disconnected exchange diagrams with brackets. These contributions correspond to the negative term in (2J. In the 
exchange diagrams we account only for s-p mixing by VJ, contributions of higher angular momenta are neglected, 
see e.g. Ref. j5j- The deformation (dashed line) is also attached only to s-p lines because it is practically saturated 
by 6s- and 6p-states. Since VJ and V x are single particle operators we evaluate each diagram solving equations for 
corresponding corrections. For example the first diagram contains on the right top leg the correction 

5^ = 2Z {mPl/2lVxlnS) \m Pl/2 ). (16) 
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To evaluate the correction we do not use a direct summation, but instead solve the equation 

(H - e)Sip x = —V x \ns) , e = e ns (17) 

for each particular \ns) state. Here H is the Dirac Hamiltonian with potential ||BJ|. Similarly the bottom left leg of 
the diagram 5 is evaluated using 

(H - e ns )8tp d = -VJ\ns) (18) 
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FIG. 5: Third order perturbation theory diagrams corresponding to Eq. 11411 . The cross denotes the T,P-odd interaction VJ 
JH3J, the dashed line denotes the lattice deformation perturbation V x W'll . and the wave line denotes the Coulomb interaction 
Vc 11311 . Multipolarity of the Coulomb interaction is shown near the line. Each diagram contributes with a coefficient shown 
before the diagram. Summation over all intermediate states \k) and \m) and over all filled states \n) is assumed. 



Diagrams with two single particle operators on the leg require a more careful treatment. For example to evaluate the 
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top right leg of the diagram 7 we first find 5ip x using Eq. I|17|l . After that we calculate Sipdx using 

(H - e ns )5ip dx = -VIH, + (ns\VJ\6iP x )\ns) . (19) 

The additional term (ns\VJ\5tp x )\ns) in the right hand side of the equation is due to the orthogonality condition 
(Sipdx\ns} = 0. The corrections Sipdx in diagrams 13 and 14 clearly do not require any additional terms in the corre- 
sponding equations. Finally the disconnected diagrams (diagrams with brackets) which originate from the second term 
in (|14fl contain the energy denominator squared. For example the bracket in the diagram 9 is equal to 2(ns\VJ\S(f> x ) 
where 

c, \ - (mpi/ 2 \Vx\ns) 

Hx = 7 72 \ m Pl/2) ■ (20) 

m \ e ns e mp l/2 ) 

To find S(f> x we calculate 5ip x at e = e ns ± S, see Eq. (|17|) . and then find S(f> x using numerical differentiation 

dip x (e ns + S)-Stp x (e ns -S) 



(21) 

8^0 



25 

As a result of the calculations we get the following T,P-odd energy correction related to the displacement x 

Ae(x) = —A— ( — ) E (n s ■ n x ) , 

as \eciB J 

A = -0.094- 0.159 + 0.080 + 0.133 + 0.141 - 0.295 - 0.257 (22) 
- 0.040 + 0.610 - 0.295 - 0.009 - 0.001 + 0.440 - 0.159 
= 0.095 . 

We remind that as is the Bohr radius, and e = |e| is the elementary charge, so the ratio of the electron EDM d e to 
ea,B is dimensionless. The results in (1221 is given is atomic units, Eq = 27.2eV. The unit vectors are ns = S/S, 
n x = x/x, where S = 7/2 is spin of Gd ion. Fourteen terms in A represent contributions of fifteen diagrams Fig. [S] 
For diagrams 13 and 14 we present only the combined contribution (+0.440). The point is that each of these two 
diagrams gives a very large contribution inversely proportional to the fine structure splitting. However these large 
contributions are canceled out in the sum of 13 and 14, so the sum remains finite even at zero fine structure splitting. 
For this reason we do not calculate 13 and 14 separately. 



IV. T,P-ODD VOLTAGE ACROSS MAGNETICALLY POLARIZED SAMPLE OF GDIG, AND 
MAGNETIZATION OF GGG IN THE EXTERNAL ELECTRIC FIELD 

In a magnetically polarized sample, according to Eq. (|22p. each Gd 3+ ion can gain energy from a small distortion 
of the lattice. The lattice has some stiffness and therefore the variation of energy per Gd ion as a function of the 
displacement x is of the form 

Ae(a:) = ]-Kx 2 - A——E Q , (23) 
2 as eaB 

where K is the effective elastic constant per Gd 3+ ion. Minimum value of Ae(x) corresponds to the shift 

as Ka B ea,B 

which is the new equilibrium position of Gd 3+ in the magnetically polarized sample. 

To find x one needs to know K. In Ref. 11] value of this constant has been estimated using the known static dielectric 
constant of GdIG e w 15 and using a rather simplistic picture of the dielectric polarization. In the present work we 
estimate K using analysis of data on infrared spectroscopy of garnets 0, |2^| . According to 0, |2^ the so called "N" 
and the "O" infrared modes are related to movement of Gd. The N mode in GdIG has energy u>n = 213cm -1 . Energy 
of the O mode is unfortunately unknown. However similar modes are known in Yttrium iron garnet: = 208.8cm -1 , 
luo = 144cm _1 . Assuming the same ratio ojjv/wo we find for GdIG: loo = 147cm _1 . Movement of Gd is not a normal 
mode of the lattice, however the splitting between ojn and ll>o is not large. Therefore to find the effective frequency 
we average between the modes 

uj = ^(u} N + wo) = 180cm _1 . (25) 
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The elastic constant is equal to 

K = (26) 
where y, is the reduced mass corresponding to movement of Gd with respect to the lattice of the garnet Gd3FesOi2, 

1 1 3 



H M Gd 5M Fe + YIMo ' 
Together with (|25|l and I12GII this gives 



/i = 77. (27) 



E 

K = 0.095^. (28) 



a 



B 



Accidentally the dimensionless coefficient in this equation is the same as that in Eq. (|22|l . The value l)28|) is by factor 
2.35 smaller than that obtained in Ref. from the dielectric constant. We believe that the present analysis is more 
reliable and gives the more accurate value of the elastic constant. Further analysis of the lattice vibrations and a new 
data on the infrared absorption would be very helpful for precise determination of the elastic constant. 

Using Eqs. (|24|l and (|28|l we obtain the value of the displacement induced by the electron EDM, x ss d e /e. 
Macroscopic electric polarization arising from the shift of Gd ions in the sample is given by P — ZexriQd, where 
nod — 1-235 x 10 22 cwr 3 is number density of Gd in GdlG. Hence, the electric field in the sample is 

E = -AttP fa l2nn Gd d e = -1.1 x 10~ 10 V/cm. (29) 



The numerical value corresponds to the current upper limit on the electron EDM For a 10cm sample this gives a 
voltage AV = 1.1 x 10~ 9 U. This value is a factor 3 greater than that obtained in |ll|. The difference is mainly due 
to the different elastic constant. 

Another effect is magnetization in the external electric field. A candidate for such experiment is gadolinium gallium 
garnet (GGG), see Ref. 0. There is no need to repeat the calculation ^l|j we need only to rescale the value. According 
to the present work the value of A in Eq. 122|) is a factor 1.27 larger than that from [TlJ and the elastic constant is a 
factor 2.35 smaller than that from 

[HI (we assume that the elastic constants for GdlG and GGG coincide). Therefore, 
the energy shift of Gd 3+ ion in GGG in the external electric field is 

Ae = (n E ■ n s ) 5.7 x 10" 22 eU , (30) 

where tie is the unit vector along the electric field, and ns is the unit vector along the ion spin. The value corresponds 
to the current upper limit on d e Q and the electric field inside the sample E — 10 kV/cm. The energy shift 
(|30|l leads to the macroscopic magnetization of the sample. The magnetization depends also on temperature and 
internal magnetic interactions in the compound. We do not discuss these points here. According to estimates the 
magnetization due to the energy shift ~ !0~ 22 eV can be measured and moreover the prospects for improvement of 
sensitivity are very good. 



V. CONCLUSIONS 



In the present work we have calculated the T,P-odd effects induced by the electron electric dipole moment (EDM) 
in gadolinium gallium garnet, GdaGasO^ (GGG), and gadolinium iron garnet GdaFesO^ (GdlG). Both GdlG and 
GGG have uncompensated electron spins on Gd 3+ ions. There are two possibilities to probe the electron EDM. 
The first one is to polarize magnetically the electron spins and to measure the induced voltage across the sample. 
According to the present calculation at the current limitation on the electron EDM, (JTJ, the induced voltage across 
a 10cm sample is 1.1 x 10~ 9 V. Another possibility is to apply an electric field to the unpolarized sample. This leads 
to the spin-dependent energy shift of each Gd ion Ae = 5.7 x 10~ 22 eU at E = lOkV/cm. This can be measured via 
macroscopic magnetization of the sample. 

In the present work we generally followed the path of 11]. However the electronic part of the pro blem has been 
considered much more accurately with account of exchange diagrams which have been neglected in |ll| . Surprisingly, 
because of accidental compensations, the final electronic effective Hamiltonian proved to be close to that from 111. 
Nevertheless the observable effects are by 3 times bigger than that in [llj . The main reason for this is in the different 
lattice elastic constant. To determine value of the constant in the present work we have used data on infrared 
absorption of garnets. A new data on the infrared absorption and further analysis of the lattice vibrations would 
be very helpful for more accurate calculations of observable effects. However the main problem which bothers us 
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is the electronic part. In the present work we have used a "jelly model" smearing eight oxygen ions surrounding 
gadolinium ion over a spherical shell. It is a sensible simplification. However we found strong compensations between 
contributions of different exchange diagrams. Clearly the compensations are related to the Schiff theorem. In this 
situation any simplification causes questions and, in our opinion, a further calculation of the effect with account of 
real 3D geometry is necessary. 
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